Glass and Ceramics 


Vol. 63, Nos. 5-6, 2006 


COMPUTERIZED MONITORING SYSTEMS 


UDC 666.1.022:658.562.012.7 

STATISTICAL CONTROL METHODS 
IN GLASS BATCH PREPARATION PROCESS 


V. V. Efremenkov 1 and K. Yu. Subbotin 1 


Translated from Steklo i Keramika, No. 5, pp. 3-5, May, 2006. 


Statistic methods for quality control of raw materials and glass batches are considered. An algorithm for statis¬ 
tic control of the batch preparation process is proposed. 


Statistical control of any technological process consists 
of correcting input and output process parameters based on 
sampling data to ensure adequate quality of finished product 
[1]. The adjustment of process parameters in the simplest 
cases can be implemented without participation of the per¬ 
sonnel, i.e., by automatically maintaining the constancy of 
one or several parameters with a certain accuracy. Such sys¬ 
tems may use deviation control, disturbance control, or a 
combined principle. 

In automated or so-called “man-machine” systems an ad¬ 
justment is performed with human participation based on sta¬ 
tistical analysis of output parameters of the process. 

The present study considers the use of statistical control 
methods in the process of glass batch preparation for sheet 
glass production at the Borskii Glass Works JSC. 

The statistical control method consists of the following 
operations [2]: 

- construction of a model of the batch preparation pro¬ 
cess; 

- identifying processes and operations; 

- analysis of the effect of the quality of raw materials 
and precision of individual operation on the batch quality; 

- analysis of the efficiency of methods used to control 
the quality of raw materials; 

- analysis of statistic control methods at production fa¬ 
cilities; 

- estimate of manageability of technological processes 
of material pretreatment and batch preparation; 

- developing proposals for statistic control and regula¬ 
tion methods; 

- experimental verification of the efficiency of propo¬ 
sals and implementation of these results in production. 


1 Stromizmeritel’ Joint-Stock Company, Nizhny Novgorod, Russia. 


The process of preparing glass batch in the general case 
consists of five interrelated operations [3] and constitutes a 
technological chain of transformations including acceptance 
and storage of raw materials, control and treatment of glass 
batch components, preparation of the batch, and its storage 
and transportation to the glass-melting furnace. 

The input control of raw materials performed by the cen¬ 
tral plant laboratory and technical services is needed to es¬ 
tablish the compliance of the quality of materials with stan¬ 
dard requirements for the following parameters: 

- chemical homogeneity and constancy of chemical 
composition; 

- content of iron oxides; 

- homogeneity and constancy of granulometric compo¬ 
sition. 

To estimate the dynamics of the fluctuations of material 
parameters from one batch to another, statistic acceptance 
control based on mathematical statistics is used. The quality 
is estimated based on the probability distribution density of 
the variations of the analyzed parameter for each material. 
One of the methods for estimating distribution density is the 
construction of a histogram in order to identify the probabi¬ 
lity distribution density law and to obtain the mathematical 
expectation and the dispersion of the controlled parameter 
values in the material. 

By way of example, let us consider histograms con¬ 
structed to estimate the quality of sodium sulfate (the content 
of the main material) in different delivery lots [4], According 
to technological standard, the content of sodium sulfate con¬ 
verted to dry matter should be within the interval: 

V max - X min , 

where Y max = 99.0% and X mm = 97.3% are the maximum and 
the minimum content of Na,S0 4 in the material. 

139 

0361-7610/06/0506-0139 © 2006 Springer Science+Business Media, Inc. 





V. V. Efremenkov and K. Yu. Subbotin 


140 





Fig. 1. Histograms of the con¬ 
tent of sodium sulfate in raw 
material with different forms 
and distribution center posi¬ 
tions: N) number of samples. 


Based on the distribution of the content of the main sub¬ 
stance on the histograms (Fig. 1) it is possible to find out 
whether the delivered material has a satisfactory quality 
level. For this purpose, based on determined tolerance, the 
number of intervals controlled, and the chosen number of 
samples (our example uses 100 samples of material), the dis¬ 
tribution form and width are determined with respect to the 
admissible tolerance width. 

The distribution form in Fig. la is satisfactory, as the left 
and the right sides are nearly symmetrical with respect to the 
distribution center, the distribution center coincides with the 
tolerance field center, and the distribution width is equal to 
the tolerance width. As the distribution domain coincides 
with the tolerance width, there is a probability of getting 
samples whose parameters are outside the established limits. 
In this case, it is advisable to tell the supplier of the need to 
stabilize the quality of material delivered. 

In Fig. lb the distribution form and center are shifted to 
the left, indicating a risk of receiving material with a de¬ 
creased content of the main substance, below the bottom to¬ 
lerance limit, which is inadmissible. 

In Fig. lc the distribution center coincides with the tole¬ 
rance field center, but the distribution width exceeds the to¬ 
lerance field width on both sides. Such material cannot be 
used in batch preparation, since it does not satisfy the estab¬ 
lished requirements on a stable content of the main sub¬ 
stance. 


TABLE 1 


Parameter 

Insoluble 

residue 

Soda 

Sulfates 

Carbo¬ 

nates 

Moisture 

Tempe¬ 

rature 

Insoluble residue 

i 

-0.40 

-0.19 

-0.54 

0.07 

0.34 

Soda 

- 

1.00 

0.22 

0.06 

0.05 

-0.25 

Sulfates 

- 

- 

1.00 

0.04 

-0.17 

-0.03 

Carbonates 

- 

- 

- 

1.00 

-0.06 

-0.30 

Moisture 

- 

- 

- 

- 

1.00 

-0.14 

Temperature 

- 

- 

- 

- 

- 

1.00 


Thus, a histogram makes it possible to recognize the 
quality of material delivered. Flowever, it does not provide 
numerical data on the distribution center, dispersion values, 
and symmetry between the left and the right distribution 
sides. To estimate the distribution center, the arithmetic mean 
and the dispersion, i.e., the mean quadratic deviation, are cal¬ 
culated. 

In the cases when significant disturbances occur in the 
process or when correctives are introduced in the operating 
mode of the equipment, point estimates, such as the arithme¬ 
tic mean and the standard deviation, are not sufficiently in¬ 
formative, since they cannot identifying the changes that 
have occurred. A stationary technological process may be 
also disturbed by other reasons: a variation in the properties 
of material used, design modifications of the technological 
equipment, etc. To identify such effects, we analyze the cha¬ 
racteristic specifics of the probability distribution density of 
analyzed parameters using nonparametric approximation of 
the corresponding statistic distribution law. 

The statistic analysis of the accuracy and stability of the 
batch preparation process, in addition to the analysis of the 
quality of material supplied, includes the solution of the fol¬ 
lowing problems: 

- comparing the mean values of the batch parameters 
and finished glass parameters to identify their random differ¬ 
ences or regular correlations; 

- comparing the dispersion (standard deviations) of the 
batch composition parameters; 

- determining the correlation coefficient between the 
batch composition parameters, 

- identifying the effect of the batch parameters on the 
variation of finished glass parameters, etc. 

In statistic analysis of glass batch quality it is essential to 
know the independence of individual parameters, which sim¬ 
plifies the control process. To reveal these correlations, we 
carried out the correlation analysis of mean daily parameters 
of the batch quality for three years of performance of the 
batch-preparing division at the Borskii Glass Works. The 
matrix of the estimated pairwise correlation coefficients 
characterizing the tightness of the relations between the 
batch parameters is shown in Table 1. 

Since the modulus of all estimated pairwise correlation 
coefficients does not exceed 0.7, the relationship between the 
batch quality parameters is considered weak [5] and can be 
neglected in developing corrective measures to stabilize the 
batch parameters. 

The statistic parameters of batch preparation quality for 
the same period are given in Table 2. 

It can be seen that the highest variability estimated based 
on the variation coefficients is observed in the content of sul¬ 
fate, batch moisture, and batch temperature. Special attention 
should be paid to stabilizing such parameters. 

At the same time, the table data do not provide an esti¬ 
mate of the variation of these parameters in time and do not 
allow taking preventive measures. Therefore, to estimate the 
accuracy and stability of the batch preparation process, plots 
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TABLE 2 


Mass content, % Batch Batch 


Parameter 

insoluble 

residue 

soda 

sulfates 

carbo¬ 

nates 

moisture, tempera- 
% ture, °C 

Mean value 

61.81 

17.88 

0.81 

20.52 

4.45 

45.80 

Standard value 

0.30 

0.26 

0.20 

0.12 

0.17 

1.55 

Variation 

0.49 

1.45 

24.70 

0.56 

3.80 

3.40 


coefficient, % 


of the variation of the quality parameters in time are con¬ 
structed, and the precision of the process is calculated from 
the formula 

K p = 6S X /AX, 

where K is the precision coefficient of the process; AX is the 
analyzed parameter tolerance; and S is the mean quadratic 
deviation. 

The deviation of the content of soda in the batch from the 
estimated values according to the results of statistical control 
in 2004 are shown in Fig. 2. It can be seen that the variation 
of soda content in time is characterized by the existence of a 
linear trend. The mathematical description of the tend in 
given in the field in Fig. 2 where the parameter “x” means the 
time in the day. On the average the content of soda during the 
year is 0.1% higher than the estimated value, which calls for 
attentive monitoring of the dosing-mixing lines. 

Similar plots are constructed for other glass batch com¬ 
ponents and parameters (temperature, moisture, content of 
insoluble residue). 

If a process is characterized by several parameters, for 
clarity the information is represented in a simple compact 
form of a circular radar plot [6]. 

The obtained statistic analysis data enable us to imple¬ 
ment statistic control of the batch preparation process. The 
statistic control algorithm can be represented as follows: 

- determination of the deviation of process parameters 
from prescribed standards using statistic methods for moni¬ 
toring the precision and stability of the processes; 

- choice of the most significant factors influencing deci¬ 
sion making using a priori data on the process; 

- evaluating these factors by correlation and multiple re¬ 
gression analysis; 



Days in 2004 

Fig. 2. Deviation of soda content in the batch from estimated values. 

- improving the process by introducing corrective mea¬ 
sures; 

- obtaining the confirmation of results and processing 
new statistical data on the technological process; 

The described statistical control methods can be effi¬ 
ciently implemented using personal computers and can be in¬ 
tegrated in the quality management system of the batch pre¬ 
paration division for producing various glass products. 
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